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QCD corrections to non-standard WZ+jet production 
with leptonic decays at the LHC 
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We discuss the impact of anomalous WW7 and WWZ couplings on WZ+jet production at next- 
to-leading order QCD, including full leptonic decays of the electroweak gauge bosons. While the 
inclusive hadronic cross sections do not exhibit any particular sensitivity to anomalous couplings once 
the residual QCD scale uncertainties are taken into account, the transverse momentum distributions 
show substantial deviations from the Standard Model, provided that the anomalous vertices are 
probed at large enough momentum transfers. 
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INTRODUCTION 



Searches for anomalous trilinear electroweak couplings 
at hadron colliders is an active field of research to dis- 
cover new interactions beyond the well-established Stan- 
dard Model of particle physics (SM), see, e.g., [l| for 
very recent D0 results. Crucial for revealing modifica- 
tions of the electroweak sector is, without any doubt, ade- 
quate experimental and theoretical precision of measure- 
ment and simulation, respectively. Since corrections from 
Quantum chromodynamics (QCD) are typically sizable 
at hadron colliders, they have to be properly included 
in phenomenological studies in order not to misinterpret 
QCD dynamics for beyond-the-SM (BSM) interactions 
when data is analyzed. Hence, an important step to- 
wards resolving these QCD-related issues is studying im- 
portant production processes at least to next-to-leading 
order (NLO) QCD precision. 

The NLO QCD corrections to anomalous W Z and 
production have already been discussed in the literature 
in detail, e.g. in Rcfs. @, 9. At the LHC, with its 
large center-of-mass energy of expected 14 TeV, addi- 
tional jet radiation from accessing gluon-induced chan- 
nels with large parton luminosities is substantially im- 
portant for these processes, which are exclusively quark- 
induced at LO. Moreover, the additional parton emission 
represents a theoretically indispensable real emission con- 
tribution to the infrared-safe NLO QCD diboson cross 
section via the Kinoshita-Lee-Nauenberg theorem Q ■ By 
probing the gluon-induced channels at small momentum 



fractions with respect to the incoming protons, the one- 
jet-contribution becomes comparable to the leading order 
(zero jet) cross section for typical selection criteria*, 



a(pp — > 3 leptons +f T + jet) 
a(pp ->■ 3 leptons + f T ) 



1 . 



(1) 



Hence, for inclusive measurements, the perturbative un- 
certainty of NLO W Z production is heavily affected by 
the one jet contribution, which is leading order (LO) in 
the strong coupling-expansion. 

Since the gluon-induced real emission contributions 
kinematically obstruct anomalous coupling measure- 
ments, one often reverts to jet vetos, at least in some 
way, when performing phenomenological studies, e.g. 
in Ref. Hence, considering exclusive NLO produc- 
tion 0, one is able to increase the sensitivity of the 
WZ cross section to anomalous couplings, which we de- 
fine to be cx(non-SM)/<7(SM) in the following. From a 
perturbative QCD point of view, however, this strategy 
does not necessarily give rise to improved QCD precision 
of total cross sections, even if improved renormalization 
and factorization scale dependencies might superficially 
suggest so. To arrive at seemingly stable NLO results, 
one subtracts a dominant and unreliable LO-a s contri- 
bution. Since this very contribution is sensitive to the 
applied hadronic cuts, one is able, given a relation anal- 
ogous to Eq. ([T]), to tune the real emission contribution 
to milden the factorization and renormalization scale de- 
pendencies at the integrated cross section level L This 
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*We refer to the considered processes pp — > l~ &i£' + £'~ + jet and 
pp — > l + vil IJr £'~ + jet, where 1,1' denote light but distinct lepton 
flavors, as W Z + jet production in this paper. 

T Indeed, this is the case for e.g. W7+jet production, where a more 
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signalizes the necessity to also take into account the dif- 
ferential QCD corrections to diboson production in asso- 
ciation with a hard jet, to add NLO precision to the veto 
performance. The quantitative dependence of the QCD 
corrections on the anomalous couplings in this context 
supplements the necessary additional information to re- 
cent theoretical leading order analysis Q. 

Eq. (TTJ also raises the question whether there is enough 
sensitivity to anomalous couplings left at larger inclusive 
rates, which could supplement the powerful "traditional" 
measurement strategies, that exploit the well-known clas- 
sical radiation zeros in the qQ —> W*y, WZ amplitudes in 
the SM @. Given Eq. (fTJ), the sensitivity to anomalous 
couplings in inclusive phenomenological analysis crucially 
depends on the one jet-inclusive cross section. A reliable 
assessment of the anomalous couplings' impact is then 
only possible if QCD corrections are properly included 
as they have turned out to be particularly sizable for the 
diboson+jet cross sections in Refs. |7Hlfj||. 

To contribute to resolving the above issues, we report 
in this paper on the calculation of non-standard WZ+jet 
production at NLO QCD. We compute the WZ+jet cross 
sections and distributions as functions of the anomalous 
couplings' parameters by means of a fully-flexible Monte 
Carlo program that has been designed for the purpose of 
this work. Full leptonic decays of the W and Z bosons 
are included throughout. The code will be publicly avail- 
able with an upcoming update of Vbfnlo [ll| , and is also 
capable of performing Tevatron calculations. We give de- 
tails on the computation in SecHH and Sec. lIIIl is devoted 
to the numerical results. In particular, we compare the 
impact of anomalous couplings to the scale uncertainty 
of the SM cross section for inclusive selection cuts and 
for a cut choice that mimics genuine WZ events. For 
the latter, we discuss the anomalous couplings' impact 
on the differential cross sections in more detail. Sec. IIVI 
closes with a summary of this work. 



II. DETAILS OF THE CALCULATION 

The strategy and the methods we apply to perform 
the NLO calculation in a fast and numerically stable way 
have already been presented in detail in Refs. jsLHol]. We 
therefore only sketch the numerical implementation to 
make this work self-consistent. To arrive at the numerical 
results of Sec. IIII1 we use a modified version of the Monte 
Carlo code of Ref. [10| , which builds upon the Vbfnlo 
framework We evaluate the LO matrix elements for 
pp — > 3 leptons + f + jet at 0(a s a 4 ) via Helas rou- 
tines [13], which are set up with MadGraph [Is} . The 
virtual and real emission matrix elements are computed 



thorough investigation reveals sophisticated cancellations of scale 
dependencies for the NLO exclusive cross sections, rendering the 
perturbative stability unreliable. 



numerically usingin-house routines along the lines de- 
scribed in Refs. [1, Il0| . Algebraic bookkeeping of the 
infrared singularities is thereby performed by applying 
the method of Catani and Seymour fli| . Throughout 
the computation, the anomalous trilinear couplings are 
included via purpose-built Helas routines, which we de- 
sign using an in-house framework that employs Feyn- 
Rules [l5[ . The anomalous trilinear vertices follow from 
the most general Lorentz-, CV-, and QED-invariant ef- 
fective operators up to dimension six [2, EH , 

Cwwj = ~ ie[wl u W^A v - WlA v W^ 

A (2) 
+ Kj WlW„F^ + -^-wLw>fF vX ] 
m w 

for the anomalous WWj vertex, and 

Cwwz = - iecote w [gf (W^WA* - W^W^) 

+ n z wlw v z^ + \w\^z vX \ 

m w 

(3) 

for the anomalous WWZ vertex. In Eq. ([3]) 9 W denotes 
the weak mixing angle. W^, A^, and Z 11 denote the W, 
the photon and the Z boson. The field strength ten- 
sors are = d^A" - d v VP , Z» v = d»Z v - d"Z», 
and W^ = d»W v - d v W», as usual. We do not in- 
clude CV- violating operators since they are already heav- 
ily constrained by measurements of the neutron's electric 
dipole moment (Ref. LT3] ) • We thus can focus on the 
W~ Z+jet production for brevity, and our findings gen- 
eralize to W + Z+ jet production in a rather straightfor- 
ward fashion. The W + Z+jet cross sections are, however, 
larger by a factor 1.5 due to testing the protons' valence 
quarks. 

We apply the gauge constraints, which were used for 
the combined analysis of LEP data in Ref. [HI , 

K z = 9\ - (« 7 - 1) tan 2 6 W , X z — A 7 . (4) 

This analysis also gives rise to the currently most strin- 
gent bounds on anomalous trilinear couplings, 

gf = 0.991±°;° 22 , « 7 = 0.984+°;°f 7 , 

A 7 = -0.016t°;{g*, (5) 

at 68% confidence level. Most constraining bounds at 
hadron colliders have only recently been obtained from 
direct measurements of the WWZ vertex at the Fermilab 
D0 0]. 

Unitarity requires the anomalous parameters 
{gf , Kz, Az, k 7 , A 7 } to be understood as low-energy 
form factors. Their precise momentum dependence is 
determined by the BSM Lagrangian. Avoiding any 
particular assumption about how the BSM interactions 
might look like (except for participation in SM inter- 
actions of some kind), we use the conventional dipole 
parametrization of the form factors, Ref. [2| . Introducing 
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new parameters, which rephrase the modifications of 
Eqs. @ and (J3J around the SM Lagrangian, 



(Afff , Ak z , Ak 7 ) = (gi,K, z ,K 7 ) - 1 



(6) 



the dipole profile is given by 



(A.9i , Ak z , X z , Ak 7 , A 7 ) 



(1 



/A 2 



(7) 

where towz denotes the invariant mass of the decay- 
ing WZ or IV 7 pair, and A is the scale where new 
physics enters the picture. The momentum dependence 
of Eq. ([7]) guarantees vanishing contributions from phase 
space regions where the naive momentum-independent 
parametrization violates unitarity. In addition, given the 
steeply falling parton luminosities for large parton mo- 
mentum fractions, a well-behaved theory yields vanish- 
ing contributions in this particular region, so that the 
parametrization of Eq. ([7]) is adequate for phenomeno- 
logical studies. In the following, we choose A = 2 TeV, 
which is also a benchmark point of various experimental 
studies (e.g. in Ref. 0, Q). 



III. NUMERICAL RESULTS 

For the numerical results, we use CTEQ6M parton dis- 
tributions pj] and the CTEQ6L1 set at LO. We choose 
m z = 91.188 GeV, m w = 80.419 GeV and G F = 
1.16639 x 10 -5 GeV -2 as electroweak input parameters 
and derive the electromagnetic coupling a and the weak 
mixing angle from Standard Model-tree level relations. 
The LO and NLO running of the strong coupling a s is 
determined by a\°{m z ) =0.130 and af°(m z ) =0.118 
with five active flavors, respectively. We fix the center- 
of-mass energy to 14 TeV for LHC collisions, and we sum 
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FIG. 1: Scale dependence of the SM cross section at LO and 
NLO for the set up described in Sec. IIIII In addition to the 
fixed-scale variation, we also plot the variation for renormal- 
ization and factorization scales chosen to be the maximum 
transverse momentum of the charged leptons. Varying fi by a 
factor two around £ = 1 reduces the uncertainties from 24% 
(23%) at LO to 9% (10%) at NLO for fj, = £ x 100 GeV 
(fi = £ maxp T ). 
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FIG. 2: Total LO cross section contours in fb for W~ Z+jet 
production for the cut and parameter choices as described in 
the text. We use a fixed-scale choice jj,R = /if = 100 GeV. 
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FIG. 3: Total NLO cross section contours in fb for W Z+jet 
production analogous to Fig. [2] 

over the light lepton flavors in the final state, which we 
treat as massless. 

Additionally, we neglect processes subject to Pauli- 
interference, i.e. we quote results for pp — > e~t? e /i + /i~j + 
X plus pp — > /j,~P tl e + e^j + X in case of W~ Z + jet 
production *. Lowering the available LHC energy to 7 
TeV yields a too low rate to be of phenomenological im- 
portance given the scheduled integrated luminosity of 1 
fb . Since anomalous couplings alter the high-energy 
phenomenology, WZ +jet production from Tevatron col- 
lisions does not yield any notable sensitivity in the al- 
lowed parameter range. 

The CKM matrix is assumed to be diagonal, which is 
an excellent approximation for LHC cross sections [To| . 
Jets are recombined via the algorithm of [20] from par- 
tons which fall in the pseudorapidity range of \r)\ < 5. 
The jet resolution parameter is chosen to be D = 0.7, 



* Since Pauli-interference is a small effect for our process, we could as 
well multiply our results for the total cross section and differential 
distributions by a global factor 1.5 to account for the processes 
with identical final state lepton flavors. 
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FIG. 4: Total K factor contours for W Z + jet production 
analogous to Fig. [2] 
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FIG. 5: Total NLO cross section contours in fb for W~ Z +jet 
production for the cut and parameter choices as described in 
the text, with the additional requirement of Eq. (|13[) . We 
again choose Ag ± ' = 0. We choose fin = /if = maxp^. 

and the jets are required to lie in the rapidity range ac- 
cessible to the detector 



\ yj \ <4.5. 
For the charged leptons, we request 
M < 2.5. 



(8) 



(9) 



We choose a very inclusive cut set-up in order to analyze 
the anomalous couplings' impact on the phenomenology 
of WZ + jet production over a wide range of accessible 
phase space at the LHC. We require 

p T > 30 GeV , p T > 25 GeV , f T > 25 GeV , (10) 

and, for the jet-lepton and lepton-lepton separation in 
the azimuthal angle-pseudorapidity plane 



Ret = ^(A^ + A^ 2 ,,) >0.2, %>0.2, (11) 

for all charged leptons and reconstructed jets. 

In Figs. [2] and [3J we representatively scan the NLO 
cross section for the fixed-scale choice [Ir = /ip = 



FIG. 6: Total K factor contours for W Z + jet production 
analogous to Fig. [5] 

100 GeV over the allowed parameter range for the anoma- 
lous couplings of Eq. [5] with the additional requirement 



z,o 



0. The total K factor contours. 
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(12) 



are depicted in Fig. 01 In particular, the anomalous pro- 
duction cross section mostly depends on the dimension 
six operators of Eqs. © and ([3]), which are not present 
in the SM Lagrangian. The QCD corrections are siz- 
able over the whole range of anomalous parameters al- 
lowed by LEP measurements. They are especially im- 
portant for values close to the SM. This follows from 
the QCD corrections, which turn out to be particularly 
large around the px thresholds and decreasingly impor- 
tant for larger transverse momentum values in the px 
distributions' tails. The corrections for our cut choices 
are relatively large for small values of p T 150 GeV 
(K ~ 1.3), where the distributions are unaffected by the 
anomalous couplings (Fig. [S]). The decreasing K factor 
when approaching the boundaries in Fig. [4] can be un- 
derstood accordingly: For large values of the anomalous 
couplings, the cross section is altered already at LO with 
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FIG. 7: Scale dependence of the SM cross section at LO and 
NLO for the inclusive set-up described in Sec. Mil with the 
additional requirement of Eq. (|13|l . Analogous to Fig. [Tj the 
scale uncertainty is reduced from 24% (23%) at LO to 10% 
(11%) at NLO for fj, = £ x 100 GeV (fj, = £ maxp^). 
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respect to the SM, so that its increase from QCD cor- 
rections is less significant. This implies that the NLO 
QCD corrections decrease the sensitivity to anomalous 
couplings. 

Confronting the non-SM cross sections with the per- 
turbative uncertainty of 0(10%) for SM-like production, 
Fig. [TJ it is apparent that for our inclusive cuts the impact 
of the anomalous couplings entirely drowns in the resid- 
ual QCD scale uncertainty at NLO. In Fig. [TJ we also 
plot the scale variation for the intrinsic scale max pfp, 
which is via Eq. ([7]), related to the characteristic scale 
of the anomalous couplings. Even for hard events with 
maxpj > 150 GeV, the impact of the anomalous pa- 
rameters is not apparent, which also explains the small 
percent-level deviations of the NLO cross section over the 
allowed parameter range in Fig. [31 The vanishing sensi- 
tivity arises from only small momentum transfers in the 
anomalous trilinear vertices due to our selection criteria: 
the jets recoil predominantly against the collinear WZ 
pair, which is an anomalous couplings-insensitive kine- 
matical configuration. A straightforward way to induce 
considerably larger momentum transfers while reducing 
the contributions from anomalous couplings-insensitive 
graphs, where both the W boson and the Z boson couple 
to the quark legs, is therefore requiring a large separation 
of the identically charged decay leptons. Note, that this 
reflects the kinematics of exclusive diboson production 
with the W and Z recoiling against each other at LO. A 
convenient choice is 



SM NLO uncertainty fi = (0.5 ... 2) maxpy 
Ak° = -0.003, \° = -0.06 
Ak° = -0.003, X° = -0.03 




100 150 200 250 300 350 400 450 
maxp f T [GeV] 



FIG. 8: NLO inclusive distribution of the leptons' maxi- 
mum transverse momentum maxp^ for various values of the 
anomalous couplings. The anomalous parameters which are 
not quoted are chosen to be zero. The solid line lies within 
the SM uncertainty band, and the plotted distributions are 
affected by scale uncertainty bands of equal width. For the 
shown distributions we have applied the cut of Eq. (|13|) . 

deviations from the SM phenomenology (Eq. ([7])) at LO. 
Hence, perturbation theory forces us to consider inclusive 
production to give reliable predictions at NLO; for our 
inclusive cut choices with respect to hadronic activity 
our numerical results are not affected by the mentioned 
pathologies and can be considered stable except for the 
residual scale uncertainty of about 10%. 



R e±e ,± > 1.5. (13) 

While the cross section decreases by approximately 20%, 
Fig. O reveals cross section deviations due to anomalous 
couplings of order 5% compared to the SM. Although 
this is still comparable to the cross section's scale de- 
pendence plotted in Fig. [Jj its increase compared to the 
SM entirely results from the large-p^ phase space region, 
Fig. [5J In this region, we find substantial deviations in 
the maxpj shape, which can be well outside the SM scale 
uncertainty for larger values of the anomalous couplings. 
Hence, provided a sufficiently large momentum transfer 
in the anomalous vertices, the WZ+jet production which 
is vetoed in pp —> WZ analysis, exhibits potential sen- 
sitivity to anomalous couplings via fits to the px distri- 
butions. However, the total cross section becomes tiny 
and the sensitivity to anomalous couplings decreases by 
including the NLO corrections as can be inferred from 
Figs, gl and M 

To this end, it is worthwhile to briefly comment on 
one-jet exclusive WZ + jet production. In particular, 
it has been shown that the exclusive diboson+jet cross 
sections are seemingly stable. In (Toj . however, it was 
also shown that applying an additional fixed-pT jet veto 
on the second reconstructed jet at NLO yields a poor 
perturbative reliability. In particular, applying such a 
veto can lead to negative bins already for modest scales 
of about 100 GeV. Indeed, this is the phase space region 
where anomalous couplings give rise to well-pronounced 



IV. SUMMARY 

We have discussed the impact of anomalous trilinear 
couplings on NLO QCD W~ Z + jet production at the 
LHC, including leptonic decays. We do not find any 
significant deviations of differential cross sections, unless 
we induce sufficiently large momentum transfers in the 
trilinear vertices. This can be realized by requiring 
back-to-back WZ pairs. The resulting modifications 
are characterized by large transverse momenta, and are 
well-outside the SM scale uncertainty that is intrinsic 
to our NLO QCD computation for large values of the 
anomalous couplings in the allowed range by LEP. 
In general we find that the differential cross sections' 
sensitivity to anomalous couplings decreases when 
including inclusive NLO corrections (Figs. 2] and [5]), 
while the QCD corrections do not exhibit any particular 
dependence on the anomalous parameters. Performing 
precise measurements in a full hadronic environment, 
for which our calculation is relevant, requires a careful 
analysis, taking into account all systematic effects, 
ranging from showering to detector effects. This is 
clearly beyond the scope of our calculation. Our Monte 
Carlo code will be publicly available with an upcoming 
update of Vbfnlo. 
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